Heteropolymers in a Solvent at an Interface 
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Exact bounds are obtained for the quenched free energy of a polymer with random hydrophobicities 
in the presence of an interface separating a polar from a non polar solvent. The polymer may be 
ideal or have steric self-interactions. The bounds allow to prove that a "neutral" random polymer is 
localized near the interface at any temperature, whereas a "non-neutral" chain is shown to undergo 
a derealization transition at a finite temperature. These results are valid for a quite general a priori 
probability distribution for both independent and correlated hydrophobic charges. As a particular 
case we consider random AB-copolymers and confirm recent numerical studies. 
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The statistical behavior of heteropolymers has been in- 
tensively studied in recent years 0-0] . They model ran- 
dom copolymers 0,0] and to some extent protein fold- 
ing |Q. For example, a chain composed by hydrophobic 
and hydrophilic (polar or charged) components in a po- 
lar (aqueous) solvent evolves toward conformations where 
the hydrophobic part are buried in order to avoid water, 
wheras the polar part is mainly exposed to the solvent 
0]. This is what commonly happens to proteins and 
makes them soluble in aqueous solutions. However other 
proteins (e.g. structural proteins) are almost insoluble 
under physiological conditions and prefer to form aggre- 
gates 0]. Many of the proteins which are insoluble in 
water are segregated into membranes which have a lipid 
bilayer structure 0]. Membrane proteins have a biologi- 
cal importance at least as great as those which are water 
soluble. Usually one distinguishes integral and non in- 
tegral membrane proteins according to whether the pro- 
tein is most immersed in the lipid bilayers or simply an- 
chored to the membrane respectively (in the latter case 
the protein is essentially water soluble) [0] . An analogous 
situation occurs for copolymers at interfaces separating 
two immiscible fluids. If the solvents are selective (i.e. 
poor for one of the two species and good for another), 
AB-copolymers are found to stabilize the interface 
Moreover, random copolymers have been claimed to be 
more effective in carrying out this re-enforcement action 



The simplest theoretical approach to the above prob- 
lems has been proposed by Garel et al. 0]. In the case 
of membrane proteins the finite layer of lipidic environ- 
ment is modeled as an infinite semi-space. Though a 
quite rough approximation, this is the simplest attempt 
in capturing the relevant features due to the competition 
of different selective effects |i~tj ]. 

We will study a lattice discretized version of their 
model. The nodes of an N links chain occupy the sites 



fi = (xn, . . . , Xid), i = 0, . . . , N of a d-dimensional hy- 
percubic lattice. A flat interface passing through the 
origin and perpendicular to the u = (1, . . . , 1) direction 
separates a polar (e.g. water), on the u ■ r > side, 
from a nonpolar (e.g. oil or air), on the u ■ r < side, 
solvent. The i-th monomer interacts with the solvent 
through its hydrophilic charge > (qi < means that 
it is hydrophobic) and contributes to the energy with a 
term — qiSgn(u ■ fi) jl^]. For simplicity we can associate 
the charge to the link between adjacent positions on the 
chain instead that to the single monomer. 

The partition function of the model for a chain W 
starting at position f is 
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z (r, {*}) = exp 1 *sgn(t 
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where /3 _1 = fcsT. If one sums over non- interacting 
(ideal) chains, then the lattice version of the model in- 
troduced in Ref. 0] is recovered. We will consider also 
the more physical case where steric interaction among 
monomers does not allow for multiple occupancy of lat- 
tice nodes, studied for a particular case in Ref. 0]. 

The free energy density of the system reads 
f(f, (3) = -limjv^oo j N lnZ(f,{q i }), where ~ denotes 
the quenched average over the distribution of the charges 
{qi}. In the following we will assume that {qi} are inde- 
pendent random variables having a Gaussian distribution 
of the form 
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More general cases will be treated at the end. In par- 
ticular considering charges not independently distributed 
is of interest for designed sequences as happens for real 
proteins. 

We will show that for a neutral chain (qo = 0) /(0, 0) < 
f(f, (3) with | f*| > N, in the large N limit, and for all (3. 
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The same holds also for q ^ if > l3 upper {\qo\, A) 
with (3 upper — > if ^ — ► 0. This implies that the chain 
is localized around the interface at any temperature if 
go = 0, and at sufficently low temperature if go 7^ 0. The 
proof is rigorous for the ideal chain, whereas for the self- 
avoiding case only a mild and well accepted hypothesis on 
the asymptotic behavior of the entropy is needed. When 
go =/= a rigorous lower bound on the free energy for both 
the ideal and the self-avoiding chain allows to determine 
a Pi 

ower (|g |, A) below which the chain is delocalized. 
We will first consider the ideal chain case and then 
explain the modifications necessary to extend the results 
to self-avoiding chains. 

Ideal chain. For clarity we derive the bounds in the 
d = 1 case. The general case does not contain any further 
difficulty Let us first consider initial positions far 

from the interface in the favorable solvent, x > N if 
go > or x < —N if qo < 0. Under these assumptions all 
chains remain in the same side, implying that sgn(xi) = 
1 (or sgn (xi) = —1 respectively) for all i. Upon averaging 
over the charge distribution, we obtain the free energy 
density of a walk in the favorable solvent: 



/* = ~ln2-|gb|. 

P 



(3) 



We give an upper bound to the free energy as fol- 
lows. Consider only chains made up of blobs of k 
steps, with k even. Bringing a blob in its globally fa- 
vored side leads to an energy contribution of the form 



Hj-th blob 



J2i=i ik(j-i 



so that 



Z(x = 0,{q t })> (C fe ) f exp<{/3^ 

j'=i 
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where Ck is the number of chains starting and ending 
in the origin and remaining in the same side. In the 
one-dimensional case it is easy to exactly determine Ck- 
1 / k \ 

It turns out Ck — -j I k I , so that, by using the 

2 + 1 V 2 / 

Stirling's formula, the asymptotic result Ck ~ 2 k k~* 
is found (in d dimensions C k ~ {2d) k k~~r' Q). The 
upper bound on the free energy is then: 



/(0,/3) < 



llnC fc 



k 

and, by using Eq. (^), we obtain: 
A/ = /(0, /?) — /*< hg (k, 0) = 



k 

i=i 



(4) 
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where the scaling function G is given by 



G{x) = 4=~e" 

V T X 



[1 — erf (a;) 
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G(x) is a positive decreasing monotonic function for pos- 
itive arguments. 

We consider separately the neutral and the q ^ 
cases. In the neutral case it turns out that the chain is 
always localized at the interface. In fact, if qo — we 
have 



Af<h (k,P) 
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It is easy to see that for any /? there exists a value k((3) 
such that ho(k,0) < for k > k{[3). For example at 
high temperature k{f3) ~ [ln(/3A)] 2 (/3A)~ 2 . This shows 
that at any temperature a neutral random chain is always 
adsorbed by the interface. 

In the non-neutral case with k = 2, one has A/ < if 



/3> A, 
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The limit lim^oo h qo (k,(3) = 0+ does not allow to de- 
duce the existence of a negative minimum in k, so that 
the previous argument, showing that the neutral chain 
is always localized, does not hold for qo ^ 0. Equation 
(g) proves localization at sufficiently low temperatures. 
For | g 1 <C A, it yields (3 upP er = v ^" 2 , and in the op- 



posite regime A <C |go|, A 
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the limit |go|/A <C 1 we can give a better estimate for 
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a larger blob size k = 2x\ (A/|go|) 
is fixed: 



— > 0, by considering 
where x > |?o|/A 
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We now look for a lower bound on / for all chain initial 
positions. If, for example, go > 0, the preferred side is the 
right one. Consider then the starting point at x = N — k 
(0 < k < N) and let gg, with E some subset of the last 
k steps of the walk (0 < \E\ < k with \E\ the number of 
elements in E), be the number of walks having the steps 
belonging to E in the unfavorable side. Upon defining 
G k = J2e 9e (G k < 2 fc ) one has 



Z(x = N — k, {q t }) = (2 N - Gk) e ^>=i qi - 



(10) 



so that, by using the inequality lnx > lna; and averaging 
over the charge distribution, one obtains 



a(A,g ,/3)l £ l -1 
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with o(A, q , (3) = exp [2A 2 /3 2 - 2f3q ] . This equation 
and its analogous in the qo < case show that there is a 
derealization temperature 



(11) 



such that, if /3 < (3i ower , f(\N — k\) > f * and the chain 
delocalizes. This argument can be extended to the d- 
dimensional case, in which 1 < Gk < (2d) k , Z^{qi) = 

(2d) N e^i qi and f* d = -j ln2d - M- 

The bounds we have proved above allow to conclude 
that there is a critical value j3 c such that for values of (3 
smaller than (3 C the chain is delocalized in the favorable 
solvent, while for larger values it is adsorbed by the inter- 
face, with the estimates Pi ower < Pc < Pupper (it is easy 
to verify that Pi ower < flupper)- The lower bound Jll| ) 
and the upper bound @, in the limit |go|/A <C 1, show 
the same behavior found by using both an Imry-Ma type 
argument (|] and variational approaches ■ 

Self-avoiding chain. All the results shown for a ran- 
dom chain can be readily generalized for a self-avoiding 
chain. Namely, a neutral chain is localized at all temper- 
atures, whereas a non-neutral chain undergoes a localiza- 
tion transition at some critical temperature (3 C . 

The derealization temperature f3i OW er can be derived 
exactly in the same way, since the division of walks into 
classes according to the number of steps made in the un- 
favorable solvent does not depend on the self-avoidance 
constraint. 

The upper bounds on the free energy, which allow to 
prove chain localization, requires instead some refine- 
ments with respect to the previous case. While the en- 
ergy term is computed in the same way as before, the 
entropy term is different. Firstly, the connective con- 
stant (k = 2d for a random walk in d dimensions) is 
different. We recall that the existence of the connective 
constant, k = Iyhin^oo In Sjv /A, for self-avoiding walks 
(SAW) has been rigorously established |ll| (Sn is the to- 
tal number of A-steps SAW starting from the same site) . 
The subleading correction of the form Sn — k n N''^ 1 
is widely agreed upon, although not rigorously proved 
|jl4j| . Secondly we introduce the notion of loop, follow- 
ing e.g. [ ft6f , and consider only walks made up of N/k 
blobs, each blob being a fc-loop, in such a way that 
different blobs can be embedded independently, as well 
as for a random chain. A A-loop is a A-steps SAW, 
starting and ending on the interface, which always re- 
mains in the same half-space, with the further condition 
£oi — xq2 < xn — Xi2 < xni — x 7v 2 Vi. It has been proved 
jl7j that the free energy density of loops is the same as 
for SAW, ki = limjv^oo \tiLn /A = k, where Ln is the 
number of A-loops. The subleading correction is usually 
assumed in the same form as for the number of SAW: 



These considerations are sufficient to generalize the 
previous results to the self-avoiding case, yielding the fol- 
lowing bounds for the critical temperature: 
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which do not depend on the assumption ( |12[ ) and is there- 
fore rigorous. Again, in the limit |i?o|/A « 1 a better 
estimate j3 uvper can be derived by using Eq. (|l2|): 
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Generic probability distribution. Up to now we have 
considered the hydrophobic charges as independently dis- 
tributed Gaussian random variables. Actually, the re- 
sults we have proved do not depend on this assumption. 
We will briefly sketch this in a few cases |ll| . 

The argument showing localization at any tempera- 
ture for a neutral chain holds true, both for random and 

\fk as fc — * 



self-avoiding chains, if 



The 



central limit theorem ensures this for independent ran- 
dom variables having a generic probability distribution 
with finite variance and null mean. In the non-neutral 
case, the existence of a derealization transition can be 
proved e.g. for a bimodal distribution. This corresponds 
to the more realistic case of two kinds of monomers, one 
hydrophilic and the other hydrophobic. We thus consider 
the generic bimodal distribution [jl9| : 

P ( qi ) = aS { qi - q+) + (1 - a) 6 (ft + g_) , (15) 

with q+,q- > 0. The probability distribution ( p"5| ) 
has three independent parameters, and fixing the av- 
erage charge qo — a (q + + q~) — Q— and the variance 
A = y/a (1 — a) (q+ + qJ) we are left with one free pa- 
rameter. It is interesting to report the derealization tem- 
perature ftiower, which provides a good estimate for the 
critical temperature in the previous cases: 
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Notice that in the limit of nearly neutral chain (\qo\ <C 
A) we get (3f ~ which is the same function of |go| 
and A as in the Gaussian case, suggesting the existence of 
a universal behavior. In the limit of nearly homogeneous 
chain (A <C |<7o| which implies a ~ or a ~ 1) instead 



bim 
ower 



ln 



diverges logarithmically 
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(12) 



a q+ 

1 — a q— 

in contrast with the Gaussian case. 

We consider now the case in which the hydropho- 
bic charges {qi} are not independent random variables, 
but are Gaussianly distributed with ql = qo Vz, qi qj — 
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q~ qj — . We assume Mff { = A 2 Vi, in anal- 

ogy with the non-correlated case, and also translational 
invariance along the chain for the correlation matrix: 
Mjj = b(\i — j\). One can prove that if long range cor- 
relations decay exponentially or even algebraically the 
neutral chain is again localized at all temperatures. In 
fact, by assuming an algebraic decay, b (r) ~ r~ v , it turns 

k s / 2 with 8 = min (77, 1) 



out that 



Only if 

correlations are so strong that they do not vanish along 
the chain (77 = 0), the chain does not localize at all tem- 
peratures. 

In the non neutral case the existence of the transition 
can be proved. For example the estimate of the dereal- 
ization temperature is 



ocorr ■ I 

Piower = min \ 



E\ 



<to\ 



(16) 



If the charges are positively correlated (M^ 1 > for 
i j), chain localization is more favored than in the non- 
correlated case, whereas if charges are anti-correlated 
(M-" 1 < for i ^ j) it is less favored. 

Asymmetric interface potentials. Finally we extend 
our demonstrations to the random AB-copolymers stud- 
ied by Sommer et al. j|. Their model corresponds to 
consider the following Hamiltonian: 

W = ■ n) + 6(- qi )9(u-n)}, (17) 



with the charges distributed according to Eq. (|T^) with 
a = 1/2 and q + = q_ (|A — 1| measures the poten- 
tial asymmetry). Such an AB-copolymer (Eq. (17)) is 



equivalent to a non-neutral chain in symmetric poten- 
tials (A = 1), a case that we have already discussed. We 
have proved the existence of a derealization transition 
for a neutral chain also in the Gaussian case. For both 
distributions, the derealization temperature shows the 
behavior (3i OW er — A ^ in the limit of nearly symmetric 
potentials (A ~ 1), in agreement with the scaling law and 
the numerical results found in ||. On the contrary, in 
the highly asymmetric cases (small and large A) different 
asymptotic behaviors for (3i ower occur fjs|] . 

To conclude, in this Letter we have proved several ex- 
act results on random heteropolymers in the presence of 
an interface. Namely, a neutral chain is localized at all 
temperatures, whereas a charged chain delocalizes at a fi- 
nite temperature. The results are quite general and hold 
for ideal and self-avoiding chains, Gaussian and bimodal 
distribution with independent and correlated charges. 
Furthermore, our lower bounds for the transition tem- 
perature confirm previous estimates. 
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